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Abstract. In terms of the fully relativistic screened Korringa-Kohn-Rostoker method 
we investigate the variation in the magnetocrystalline anisotropy energy (MAE) of 
hexagonal close-packed cobalt with the addition of platinum impurities. In particular, 
we perform calculations on a bulk cobalt system in which one of the atomic layers 
contains a fractional, substitutional platinum impurity. Our calculations show that at 
small concentrations of platinum the MAE is reduced, while at larger concentrations 
the MAE is enhanced. This change of the MAE can be attributed to an interplay 
between on-site Pt MAE contributions and induced MAE contributions on the Co 
sites. The latter ones are subject to pronounced, long-ranged Friedel-oscillations that 
can lead to significant size effects in the experimental determination of the MAE of 
nano-sized samples. 
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1. Introduction 

Cobalt alloys, such as CoPt or CoPd, are ubiquitous in the field of magnetic record- 
ing and of particular interest to the field of ultrafast magneto-optics pp. In terms of 
magnetic recording, increasing areal densities require decreased grain size, which in 
turn requires increasing values of magnetocrystalline anisotropy energy (MAE) to en- 
sure thermal stability of written information [2] . Currently this is achieved using CoPt 
alloys with perpendicular anisotropy. Consequently an understanding of the origin of 
the MAE in CoPt is an important practical problem. Since the magnetic properties of 
these alloys are highly sensitive to the amount and the spatial distribution of the Pt 
content, understanding the effects of alloying is an important issue. The effects on the 
magnetic properties of CoPt as functions of the platinum content have been studied 
extensively, both theoretically J3] and experimentally [U Ej. Moreover, in recent ex- 
perimental work [6] it was demonstrated that the magnetocrystalline anisotropy energy 
(MAE) of cobalt can be tuned by letting platinum impurities migrate into the cobalt 
system. Generally it is agreed that the addition of platinum to a magnetic material, 
such as Fe or Co, influences the magnetic properties, in particular, the MAE of the 
material primarily through the strong spin-orbit coupling of Pt [7j . 

The aim of the present work is to elucidate from first principles the effect on the MAE 
of bulk hep Co by the addition of platinum. To this end, we use the fully relativistic 
screened Korringa-Kohn-Rostoker (SKKR) method as combined with the coherent- 
potential approximation (CPA), which is well suited to describing substitutional alloys 
[8]. Our model focuses on Pt alloying in a (0001) atomic plane of a hep Co bulk 
system, from the case of an impurity to the case of a complete filling of the layer by 
Pt. After briefly discussing the computational methods we present the calculated MAE 
as a function of the Pt concentration and analyze the results in terms of layer- and 
species-resolved contributions to the MAE. We note that recording media are complex 
alloy systems, often containing Cr to promote grain boundary separation. It is often 
found that the maximum MAE as a function of Pt concentration is limited by, for 
example, the formation of new phases [5] or the presence of stacking faults [TUl [TT] [T2] . 
Here we are concerned only with the intrinsic enhancement of the MAE introduced by 
the Pt impurities. Remarkably, this analysis highlights the role of long-ranged Friedel 
oscillations in forming the MAE of the system. Specifically, we demonstrate a layer 
dependence of the valence charge which makes the effect of the Pt impurities long- 
ranged. This might have significant impact on the determination of the MAE of thin 
film samples corresponding to the systems studied in this work. In particular it might be 
expected to give rise to finite size effects in the MAE of granular thin films for magnetic 
recording which would become more significant as the grain size is reduced. 
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2. Computational details 

The central feature of the SKKR method is the evaluation of the electronic Green's 
function of a layered system. Here, a layered system refers to a system exhibiting two- 
dimensional translational symmetry in each (infinite) atomic plane, but in which there 
are no symmetry requirements along the third axis. From the Green's function one 
can then determine a number of physical quantities of interest, such as site-projected 
charges, spin- and orbital moments and the total energy of the system. As the method is 
well documented elsewhere [131 EH] , here we present only some details of our calculations. 
The calculations were performed within the local spin-density approximation (LSDA) 
of density-functional theory (DFT) as parametrised by Vosko et al. [15] The effective 
potentials and fields were treated in the framework of the atomic sphere approximation 
(ASA). The substitutional Pt alloying was treated within the coherent potential ap- 
proximation (CPA) [8], EE]. As the LSDA fails in predicting the orbital moment and the 
MAE for hep Co correctly, we employed a heuristic extension of the relativistic electron 
theory by the orbital polarisation (OP) correction [T71 [T8| [19], as implemented within 
the KKR method by Ebert and Battocletti [20] • The corresponding Kohn-Sham-Dirac 
equations were solved using a spherical wave expansion up to an angular momentum 
number of i = 3, although it should be noted that the OP correction was applied only 
for the i = 2 orbitals. 

The magnetocrystalline anisotropy energy was evaluated within the magnetic force 
theorem [21], in which the total energy of the system can be replaced by the single- 
particle (band) energy. Moreover, we employed the torque method [22J, making use of 
the fact that, for a uniaxial system, the MAE, K, can be calculated up to second order 
in spin-orbit coupling as 



where, in the case of hep geometry, 9 denotes the angle of the spin-polarisation with 
respect to the (0001) direction, i.e., the direction perpendicular to the hexagonal planes. 
Note that the z-axis of the (global) frame of reference in our calculations is defined to 
be parallel to the (0001) direction. Within the KKR formalism, K can be decomposed 
into site- and species-resolved contributions, 



where sf denotes the concentration of species a at site % and D\ denotes the 
corresponding derivative of the band energy. Using Lloyd's formula [23], can be 
calculated as [21] 



K = E{6 = 90°) - E(6 = 0°) 



dE 



6»=45° 



(1) 




(2) 
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where is the Fermi energy and, in case species a occupies site i, (e) and r^' n \e) 
stand for the angular momentum matrices of the single-site t operator and the site- 
diagonal scattering path operator, respectively. All these quantities are calculated at 
the direction of the magnetisation n = f 0, , corresponding to 9 = 45° in Eq. ( 1 ). 
The derivative of the t-matrix is evaluated as described in [25]. The energy integral in 
Eq. ^ can be accurately performed by sampling 20 energy points on an asymmetric 
mesh along a semi-circle contour in the upper complex semi-plane. In order to achieve 
an accuracy within 5 % for the MAE, a sufficiently dense mesh in the two-dimensional 
Brillouin zone (2D-BZ) was used to evaluate r\° l ' a \e): at the energy point closest to 
the Fermi energy, we used 5764 /c-points in the irreducible wedge of the 2D-BZ, corre- 
sponding to more than 34 000 fc-points in the full 2D-BZ. Due to the two-dimensional 
translational symmetry of the system, the MAE should be related to a 2D unit cell, 
therefore, in the following the index i in Eq. (|2| is used to label atomic layers. 

The SKKR method as applied to layered systems requires the system to be divided 
into a middle region wedged between two semi-infinite bulk regions. Adhering to this 
requirement, the effect of platinum alloying in a single atomic layer of bulk hep Co was 
investigated by considering a layered system as shown in Fig. [T] Each atomic layer 
in the semi-infinite bulk regions corresponds to pure hep Co bulk. Since the middle 
region needs to contain an integer number of unit cells and since each unit cell spans 
two atomic layers, this region consists of 2Nl hexagonal Co layers stacked along the 
(0001) direction. In one of the two central layers of the middle region, namely, in the 
one indexed by in Fig. [TJ a fraction s of the Co atoms are replaced by Pt atoms. From 
here on, this layer will be referred to as the impurity layer. It should also be mentioned 
that in this work no attempts are made to trace any structural relaxation effects of the 
hep Co lattice caused by Pt impurities. 

To take into account relaxation of the effective potentials and fields, we performed 
self-consistent calculations with Nl = 14, i.e. for 28 layers in total. One important 
consequence of the geometrical construction shown in Fig. [I] is that the calculation of K 
in Eq. ^ is confined to layers within the middle region, i.e., for — (Nl — 1) < i < Nl- 
This means that the long-ranged Friedel oscillations that arise due to the presence of 
Pt impurities are necessarily truncated. In order to safeguard against any numerical 
artefacts caused by this truncation, we increased the number of atomic layers in the 
middle region until the layer-resolved MAE converged to within about 1 % accuracy to 
the bulk Co MAE at the outer edges of the middle region. According to our calculations 
(see below), this condition requires Nl = 40, i.e., 80 atomic layers in total. We 
performed these calculations of the MAE by appending the perfect bulk potential of 
hep Co to the layers —39 < i < — 14 and 15 < i < 40, i.e., neglecting self-consistency 
effects for these atomic layers. To check the accuracy of this approach, we compared 
Z}j Co ^ for atomic layer no. 14 (with relaxed self-consistent potential) with that for atomic 
layer no. —14 (with appended Co bulk potential) and obtained that the two values agree 
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Figure 1. Sketch of the geometry of the system containing Nl hep unit 
cells, i.e., 2Nl atomic layers wedged between two perfect semi-infinite bulk Co 
systems. In the zero-indexed layer (black circles) a random substitutional alloy 
with platinum, Coi_ s Pt s , is considered. Note that the z-axis is defined to be 
parallel to the (0001) direction of the hep crystal. 



3. Results and Discussion 

To test our computational method, we first determined the MAE of bulk hep Co. Ex- 
cluding the OP correction we obtained an easy-plane magnetisation and a MAE of 6.7 
/ieV/Co atom, while including the OP correction we instead obtained an easy axis per- 
pendicular to the hexagonal Co planes and a MAE of 84.4 /ieV/Co. The latter result 
is in good agreement with the experimental value of 65.5 /zeV [26] and with the exper- 
imental easy axis being along the (0001) direction. Our result also compares well with 
that of Trygg et al. [27], who calculated K = 110 /zeV for hep Co using a full-potential 
LMTO method including OP correction. 

As described in Section 2, we performed calculations of the MAE of a bulk Co system in 
which a single layer has been substitutionary alloyed by Pt in a fraction of < s < 1. 
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The layer-resolved Co contributions to the MAE, Dl Co \ see Eqs. (2) and (3), are shown 
in Fig. [2] for s = 0.01 and s = 0.02. Remarkably, even such small amounts of Pt induce 
large fluctuations in D^ 00 ^: in the impurity layer (i = 0) and in the Co layers near the 
impurity layer (1 <| i |< 4) the relative changes of D^ ^ with respect to the bulk Co 
MAE reach 10 %. In particular, the Co contribution from layers i = ±1 is enhanced 
to nearly 94 fieV, while that those from layers i = ±2 are reduced to nearly 76 /ieV 
for s = 0.02. For layers further away from the impurity layer (| i |> 5), oscillations in 
Df ^ with rapidly decreasing amplitude can be seen. Reassuringly, the layer-resolved 
Co contributions approach the bulk Co MAE towards the outer edges of the middle 
region chosen in our calculations (i — > Nz = 40 and i — » —(Nl — 1) = —39). The mirror 
symmetry around the impurity layer, D^ Co ^ = D^°\ is also fulfilled with a high accuracy. 
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Figure 2. (Color online) Calculated layer-resolved Co contributions, to the 

MAE, see Eq. (|]), across the system shown in Fig. [TJfor s = 0.01 (red +) and s = 0.02 
(black •). The MAE of bulk Co is indicated by the solid black horizontal line. Solid 
lines connecting the symbols serve as guides for the eye. 



Our earlier studies of the MAE of impurities [281 [29], justified that the MAE is extremely 
sensitive to the presence of Friedel oscillations in the charge density. It is, therefore, 
tempting to relate Fig. [2] to the change in the valence charge on the Co atoms with 
respect to the distance from the impurity layer. From Fig. [3]^a) we can see that the Co 
atoms in the impurity layer and, in particular, the Co atoms in layers adjacent the im- 
purity layer gain some extra charge, while the charge transfer to more distant Co layers 
drops rapidly. The energy shift of the layer-resolved Co valence band position is well 
described by the layer-resolved change in the Madelung potential. (Note that within 
the ASA the Madelung potential in each atomic plane is a constant.) As is obvious 
from Fig. [3tb), an enhanced (reduced) charge at the Co sites is accompanied with a 
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downward (upward) shift of the valence states. Comparing with Fig. |2j this shift of the 
valence states correlates directly with the MAE contributions of the Co layers adjacent 
the impurity layer, but, clearly enough, the changes in the MAE contributions from 
more distant Co layers are also subject to fine details of the valence states influenced 
by the Pt alloying. 
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Figure 3. (a) Calculated valence charge on the Co atoms, Q, and (b) relative shift of 
the Madelung potentials with respect to the bulk case, AVMad, for layers — 13 < i < 13 
and for a Pt concentration, s = 0.02. Solid lines serve as guides for the eye. 



The effect on the species-resolved MAE by alloying with Pt is demonstrated for the 
whole range of s in Fig. 4, showing Df" ^ for layers < i < 4, together with the direct 



contribution of Pt, D Q 

JCo) 



(Pt) 



. In the impurity layer i = 0, see Fig. |4|(a), the Co contribution 
D^ 0> is reduced by the addition of Pt for concentrations up to about s = 0.15 and 
then enhanced for concentrations 0.15 < s < 0.50. For concentrations s > 0.50, D^°^ is 



Effect of Pt impurities on the magnetocrystalline anisotropy of hep Co: a first-principles study8 

again reduced with increasing s and at s ~ 0.85 D^ ^ even becomes negative. Note that 
_D Co ' ) for s — > 1 (not calculated here) would correspond to the contribution of a single 

Co atom in a pure Pt layer which, in general, would differ from zero. The on-site plat- 

(Pt) 

inum contribution, Dq , approaches the very small value of 0.01 meV as s — > 0, rapidly 

(pt) 

increases up to 0.30 meV at s ~ 0.5 and then saturates at 1 m 0.35 meV for larger s. 



It can be inferred from Fig. f4[b), that Pt alloying most dramatically influences the Co 
contribution at the layers adjacent to the impurity layer: D\ increases almost linearly 
from the bulk MAE at s = to about 0.7 meV at s — 1. As already seen in Fig. U the 
Co contributions _D- Co ' 1 from layers further out (2 < i < 4) decrease with increasing s 
and even becomes negative at s ~ 0.25 for i = 2 and 3. For s > 0.5 these contributions 
show a modest increase, but Df o) still remains negative. 

It is worth investigating the change of valence states projected onto the Co atoms in 
layer 1. From Fig. [5]^a) it is obvious that the valence charge at this Co atom increases 
almost linearly with s from 9.00 e to 9.16 e. This increase in the valence charge is 
necessarily accompanied by a downshift of the corresponding valence states, as charac- 
terized by the change in the Madelung potential, which is also linear s, see Fig. [5| 



The large enhancement of the MAE contribution from layer 1, D^ ^, can therefore be 
related directly to the monotonic shift of the corresponding valence states. 

While the species- and layer-resolved contributions to the MAE are very illuminating for 
a microscopic description of the variations in the MAE, from an experimental point of 
view only the MAE of the whole system can be accessed. Here, this means considering 
the MAE of the entire middle region illustrated in Fig. [T] for Nl = 40. In order to 
extract the change in this MAE induced by the Pt impurities, we define the excess 
MAE, AK(s), by subtracting the MAE of the 'unperturbed' cobalt bulk layers, 

40 

AK(s) = sD ( Pt) + (1 - 8 )Df o) + 2 D f° ] ~ % 1K Co , (4) 

i=i 

where K Co is the calculated MAE of hep bulk Co (84.4 /xeV). Note that we have 
taken into account the off-centre positioning of the impurity layer by doubling the Co 
contributions D^ ^ for i G [1,40], thus, in total, a system of 81 layers is considered. 
AK is shown as a function of s in Fig. [6j demonstrating that for small concentrations of 
platinum (s < 0.24) the addition of platinum to bulk cobalt actually reduces the total 
MAE of the system by about 80 /xeV. This is in strong contrast to the on-site contribu- 
tion of Pt, Dq P1: \ being positive for all values of s as seen in Fig. |4^a). The reduction 
in MAE for low s, therefore, stems from the decrease in the cobalt contributions D^ ^ 
with increasing s, in particular, for i = 0, 2, 3 and 4, see Fig. |4| AK becomes positive 
for s > 0.24 as the increasing on-site contribution, Dq\ gets larger weight (note that 

it is multiplied by s) and due to the large enhancement of D^ ^. At s = 1, AK = 1.4 

f Pt) 

meV, which is approximately four times the on-site platinum contribution, Dq 0.35 
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Figure 4. (Color online) Calculated species-resolved MAE contributions (a) for the 
impurity layer: D 1 ^ 1 ^ (black •) and D^ ^ 1 (blue k) and (b) for Co layers: D 1 ^ ^ (black 
• ), D^°^ (red +), D^ C °^ (green x) and D^°^ (purple ■). Solid lines connecting the 
symbols serve as guides for the eye. 



meV for s = 1. 



For nano-sized systems, it might be of interest to consider the change in the MAE per 
platinum atom in the system, K Pt , defined by 

K F M = , (5, 

and also, the change in the MAE per platinum atom added to the system, Kp t , obviously 
given by 

d(AK(s)) 

K Pt (s) = . (6) 
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Figure 5. (a) Calculated valence charge on the Co atom, Q, and (b) relative shift of 
the Madelung potential with respect to the bulk case, AVMad, for layer 1 as a function 
of the Pt concentration, s. Solid lines serve as guides for the eye. 

We obtained Kp t (s) by fitting a fourth-order polynomial to the function AK(s) in Fig. [6] 
and then finding the derivative of this function analytically. As apparent from Fig. [7| 
both Kp t and Kp t are monotonically increasing with increasing s, starting with the 
same value of about —1 meV at s = (see later). It follows directly from Fig. |6j that 
Kp t crosses zero at s ~ 0.24, while Kp t crosses zero at s ~ 0.11 (i.e. where the function 
AK(s) reaches its minimum). For a complete platinum layer immersed in bulk cobalt, 
i.e. for s — 1, Kp t = AK w 1.4 meV. A comparison with Fig. kj shows that about 
25 % of this value arises from the direct contribution of Pt, D[ P , and the rest from 
the induced contributions at the Co atoms. Interestingly, the change in the MAE by 
addition of a Pt atom to the system, Kp t , exhibits a surprisingly large value of about 
2.5 meV at s = 1. From Fig. [41(a) it can be inferred that the on-site Pt contribution 
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Figure 6. The change in the MAE, AK, of a system of 81 atomic layers, see Eq. Q, 
as a function of the Pt concentration, s. The solid line connecting the symbols serves 
as a guide for the eye. 
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Dq has nearly zero slope in this region of s, thus, this large value of Kp t stems mainly 
from an increase in D^ 00 ^ for 1 <| i |< 4 near s — 1. 
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Figure 7. (Color online) Black circles: calculated change in the MAE per Pt atom, 
Kp t , see Eq. as a function of the Pt concentration, s. The black solid line 
connecting the symbols serves as a guide for the eye. Red solid line: the change 
in the MAE per Pt atom added, Kp t , see Eq. (JsJ) , as calculated from a polynomial fit 
of AK{s) in Fig. [6] 
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In the limit s — > 0, corresponding to the case of a single Pt impurity in bulk Co, Kp t 
and Kp t should be identical, since for small s the function AK(s) exhibits, in principle, 
a linear dependence. This is fairly well confirmed by our calculations. Kp t (0) can then 
be expressed as 



K Pt (0) = D ( Pt) (0)-K l 



Co 



40 



dDf'Hs) 



-10 



ds 



(7) 



The physical meaning of the above equation is that adding a Pt impurity to bulk Co has 
two effects on the MAE of the system: the first two terms, Dq (0) — Kc , represent the 
direct contribution of a Co atom being replaced by a Pt atom, whereas the the last term 
of Eq. quantifies the induced change in the MAE contributions from the Co atoms 
that are not being replaced by Pt. Since the direct contribution is about -0.07 meV, 
see also Fig. |4|a), the value of Kp t (0) = — 1 meV can again only be explained by the 
induced Co contributions. 
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Figure 8. Calculated derivatives of the layer-resolved Co contributions to the 
MAE, Eq. Q, for the Pt concentration s = 0.02. The solid line serves as a 
guide for the eye. 



In Fig. [8] we show the approximate layer-resolved derivatives calculated as 



dDf o) (s) 



ds 



Df°\s 



3 + 1, 



D 



(Co), 



(8) 



for Sj = 0.02 with j indexing the Pt concentrations in ascending order. Clearly, this 
figure is closely related to Figs. [2] and [i) for small s, the Co contributions {Df" ^} show 



an increasing tendency with increasing s for 



1 and I i \> 7, while they decrease 
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for 2 <|i |< 5. Apparently, this latter effect overcomes the former one, leading to the 
relatively large value of K Pt (0) = — 1 meV. 

We also investigated possible effects of the Friedel oscillations on the MAE by truncating 
the sum in Eq. Q and considering the variation in Kp t against the number N of Co 
planes included in the sum, 



for 1 < JV < 40. The function K Pt {N) for s = 0.01, 0.05, 0.10 and 0.20 is shown in 
Fig. [9] For all cases, the maximum of Kp t (N) occurs at 2N +1 = 3 planes, i.e. includ- 
ing only one Co layer on each side of the impurity layer. This is because the induced 
effect on D^°* by the addition of platinum is strongly positive for all s, see Figs. [2] 
and[4^b). There is a significant minimum in the calculated Kp t (N) at 2N + 1 « 11 
planes. Comparing with Fig. [2j it is obvious that this minimum is due to the reduction 
of Di°°^ for 3 <| i |< 5 caused by the addition of Pt. Kp t (N) then exhibits a local 
maximum at 2N + 1 « 31, mostly due to the Co contributions in layers 8 <| % |< 15 
counterbalancing the Co contributions of opposite sign in layers 3 <| i |< 5. Concerning 
the overall accuracy of the calculated MAE, the effects of the Friedel oscillations remain 
significant for about 2A r + 1 < 70 planes, i.e. for 35 layers on either side of the impurity 
layer. In general, the variation in Kp t with N spans more than 1.5 meV for all values 
of s. This of course has significant implications for measuring the MAE in thin film 
samples, as up to about 2N + 1 < 40 planes, i.e. for film thicknesses d < 8 nm, the 
change in the MAE induced by the Pt impurities located at the centre of the sample, is 
expected to be extremely sensitive on the film thickness. Note that this finite-size effect 
is superimposed on and, most likely, amplified by quantum interferences arising from 
the boundaries of the finite film sample. 

It should be noted that, being a mean- field approach, the CPA completely neglects 
both structural and electronic short-range order effects. Such short-range order effects 
are likely to be most strongly pronounced for small Pt concentrations s. Therefore, our 
results in the low-s limit should be tested against another method. By employing a fully 
relativistic real-space embedded cluster Green's function technique as combined with 
the SKKR method [30J we made an attempt to test the effect of electronic relaxations 
around a Pt impurity placed in hep bulk Co. We performed self-consistent calculations 
for a cluster containing a Pt impurity and the neighboring Co atoms up to three nearest 
neighbor (NN) distances of the hep lattice (cl^n)- Note that our embedded cluster 
included 158 Co atoms around the central Pt atom, sorted out geometrically as follows: 
(i) 36, 2 x 30, 2 x 19 and 2 x 12 Co atoms in layers 0, ±1, ±2 and ±3, and (ii) 12, 56 
and 158 Co atoms within the spheres centered around the Pt site having the radii of a aw, 
2a,NN and 3<2aw, respectively. The MAE for this cluster was again calculated using the 
magnetic force theorem. By summing up all the site-resolved contributions of the MAE 




sD { Pt) + (1 - 8 )Df o) + 2 Df° ] - (2N + l)K Co , (9) 



i=l 



N 
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Figure 9. (Color online) Calculated excess MAE per platinum atom, Kp t , as a 
function of the number of planes N included in the sum in Eq. ([£} for s = 0.01 (red 
+), s = 0.05 (green x), s = 0.10 (blue ▲) and s = 0.20 (black •). 

in the cluster within the spheres as mentioned above, we obtained the values for Kp t , 
0.67 meV, -0.31 meV and -0.38 meV, respectively. Obviously, a direct comparison of 
these values with those for the layered system is hardly possible, since, according to the 
geometrical classification (i), the Kp t for the cluster refers to an incomplete summation 
over sites in the respective layers. Nevertheless, comparing with the values in Fig. [9] 
related to 2N + 1 = 3, 5 and 7 and for the concentration of s = 0.01 being closest the 
case of an impurity, we can conclude that both the trend and the magnitude of Kp t are 
in satisfactory agreement between the CPA and the real-space calculations. 

4. Summary and Conclusions 

Using the fully relativistic screened Korringa-Kohn-Rostoker method combined with 
the coherent potential approximation, we have studied the MAE of a bulk hep Co 
system in which one of the (0001) atomic plane has been alloyed with Pt for the whole 
concentration range, < s < 1. We conclude that low concentrations of platinum 
reduce the overall MAE of this system and that the origin of this reduction are the 
induced changes in the MAE contributions from the Co atoms. In the limit s — > (bulk 
Co), the change in MAE per platinum atom added is approximately —1 meV. At larger 
concentrations, the direct MAE contribution of the platinum, which is positive for all Pt 
concentrations, starts to increase, but the overall change in MAE due to the addition of 
platinum is still dominated by induced Co contributions. Interestingly, in the limit of a 
completely filled Pt layer, addition of one Pt atom increases the MAE of the system by 
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about 2.5 meV. We also investigated the effect of long-ranged Friedel oscillations and 
established a large sensitivity of the MAE on the number of Co layers included in the 
calculations. This might have a significant impact on the experimental determination 
of the MAE in thin film samples of this type of system. 
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